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Abstract 

An infinite number of conserved quantities in the field dynamics (f)t = L U{(f>) + p for a finear 
Hermitian (or anti-Hermitian) operator L, an arbitrary function U and a given source p are 
presented. These integrals of motion are the multipole moments of the potential created by <j) in 
the far-field. In the singular limit of a bistable scalar field (j) — (i.e. Ising limit) this theory 
describes a dissipative boundary motion (such as Stefan or Saffman- Taylor problem that is the 
continuous limit of the DLA-fractal growth) and can be exactly integrable. These conserved 
quantities are the polynomial conservation laws attributed to the integrability. The criterion 
for integrability is the uniqueness of the inverse potential problem's solution. 
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Motivation: Two-dimensional (2D) Laplacian growth (known also as 
Hele-Shaw, Saffman- Taylor, or the quasistatic Stefan problem) has many 
attributes of an exactly integrable nonlinear problem [1]. At the same time 
it is dissipative, which is not usually the case for exactly integrable the- 
ories. These problems of moving boundaries belong to the more general 
class of pattern formation phenomena which have been studied intensively 
for the last decade [2]. From different first-principles models describing 
pattern formation we will choose and discuss here the so-called phase-field 
model [3] derived from the basics of the thermodynamics and describing 
the kinetics of various phase transitions, such as different kinds of sohd- 
ification (from the supercooled liquid, from the supersaturated solution, 
from the electrolyte), phase transformations in liquid crystals and in other 
condensed matter systems, etc. It is known that this phase-field model in 
the appropriate singular limit is equivalent to the symmetrical Stefan prob- 
lem [4]. The following elementary (but unexpected) observation motivated 
this work: a quite general field dynamics (Eq.(l) below) appears to have 
an infinite number of conservation laws such that in the singular limit the 
appropriate interface dynamics can be an exactly integrable (dissipative) 
theory. Besides its purely fundamental interest, this result can help one 
to study of cases with nonzero (and even anisotropic) surface tension. It 
turns out that the conservation laws in the presence of surface tension are 
not lost, but are only corrected (Eq.(35) below). Thus in the framework of 
exact integrability one can eliminate the short-wave instability, which is of 
great physical and practical importance, because without this stabilization 
the problem of Laplacian growth is ill-posed [2]. 

Conserved integrals: Let us consider the following dynamics of a field 
4>{x,t) on the D-dimensional manifold M [xeM): 



where 0t means partial derivative of the (f) with respect to time t, L is a 
linear Hermitian or anti-Hermitian operator, U is an arbitrary (generally 
nonlinear) function of (p (and may also be a function of 0t and of V0), and 
p{x, t) is a given source term. 

Let us consider solutions, ipnix), of the equation 



<f>t^LU{<P)+p 








(2) 



such that the integrals 




(3) 
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are finite. 

Here if D > 1, the zero eigenvalue of L is infinitely degenerate, where 
n can be discrete or continuous. In general, n stands for the appropriate 
set of numbers designating the eigenfunctions of zero-eigenvalue of L. For 
example, if L is the Laplacian and D — 2 we have 

V'n = (4) 

Here n = 0,l,2,...,oo; and 2; is a complex coordinate on the plane 
z — X + iy. For D = 3, 

tPt^ =P;^!ll{e)e'^-/r^ (5) 



or 



^ p("'){0)e'"^r'' (6) 

Here, r, 6, and u are the polar coordinates and Pi^\d) are the associated 
Legendre polynomials with \m\ < n. Also n > in Eq.(5). 

One can immediately see that the dynamics of the integrals c„ is known 
and given by 

dCn/ dt = Ipn 0t 
JM 

= / il:nLUi(i))d^x+ I ilJnpd^x= I i/j^pd^x . (7) 

JM JM JM 

The first integral vanishes because L is Hermitian (or anti-Hermitian) : 

/ ipnLU{(l))d^x = ± [ U{(t))Lil:^d'^x = Q . (8) 

JM JM 

There are two important particular cases: 

A) p is constant in time. In this case the integrals c„'s are linear in time. 

B) p = 0. In this case the c„'s are constants in time. 

Harniltonian and dissipative processes: The dynamics (1) is quite gen- 
eral. Many important physical processes arc described by this equation 
with U = 6H/ 6(j), where (in the language of the theory of phase transitions) 
(f), and U are the free energy functional, the order parameter, and 
the external field respectively. We do not specify the nature of (f) yet: it 
can be a scalar, vector, tensor, etc. For the p — case, the corresponding 
dynamical equation is: 

* = Lf (9) 
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and can describe both Hamiltonian and dissipative processes. In the former 
case if is a Hamiltonian, in the latter one H can be a Lyapunov functional, 
i.e. a functional decreasing in time. 

An example of a Hamiltonian process is: 

0, = V X ^ (10) 

Here is a vector and L — curl. In this case the zero eigenvalue functions 
are t/j — grad W with arbitrary W, and it is easy to see that H is constant 
in time: 

dH rSH ^ ,^ f6H 6H^ ,^ 

An example of a dissipative process is: 

0, = (aV^ -b)— (12) 

Here is a scalar and L = aV^ — b, where a and b are positive real numbers. 
It is easy to sec that if there is no source (p = 0), H decreases in time: 

dH f5H ,r. rSH , SH 



dt 



f oH ^ ,n f dH . dH ,r, 

= -/{a(V^f + d^xO (13) 

MathemMical interpretation of the quantities Cn 's: For the particular case 
where Eq.(l) describes a Hamiltonian process, L is skew-symmetric, and a 
Poisson bracket { , } can be introduced. The integrals Casimirs [5], 

i.e. such functionals, that commute with all others: {c„, F} — for all F. 
Since we are mainly concerned with dissipative processes and L in this case, 
as a rule, is not skew-symmetric, we do not have such an interpretation of 
the c„'s. Wc can only say that the set of c^'s determines the projection of 
the field onto the null space of the operator L, (i.e. the space spanned 
by the cigcnfunctions ipn of L with zero eigenvalue). The knowledge of 
the c„'s provides only partial information about 0. In the source-free case, 
the dynamics of occurs outside the null space while the projection of 
onto the null-space of L is constant in time. So a great deal of information 
concerning the field is beyond our knowledge if the c„'s are all that we 
know. 

Physical interpretation of the Cn 's: If L is a second order elliptic operator, 
the quantities c„'s are the coefficients of the multipole expansion of the 
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potential ^{t,x) created by a distribution of "matter" with density (f)' — 
(f) — hm|2;|^oo (f) in the far-field: 

L^^cf)' (14) 
To see this, we write the solution of Eq.(14) as 

$(x)= j G{x,y)(t>\y)d''y (15) 

where 

LG{x,y)^5{x-y) (16) 
and expand the Green's function G{x,y) as 

G{x,y)^Y.^'^i^)My) (17) 

n 

where ijj'^ and ipn are the two linearly independent solutions of the equation 
Li/j = corresponding to a given n. Assuming the commutativity of the 
sum and the integral, we finally get 

^t,x)^J2<^n{t)i^'n{x) (18) 

n 

Let us take the Laplacian cases L = for and i?"^ and assume the 
field 4>'{x) to be localized so that it is zero in the far- field. In the 2D case 

m,z)^Re^Cn{t)z-'' (19) 

n 

where 

Cn(t) = / z''(l){t;x,y)dxdy, (20) 

and in the 3D case 

'^it,x) = Y.ll^nm{t)Pt\0)e'"^/r^^^^^ (21) 

n m 

here x stands for a vector-position and 

Cnm{t) = / P^^^ {9) e'^^r^ cj>{t; r, 9, u) d'x (22) 

Now the question arises: how much can we say about the density of 
"matter" 0(t, x) from the knowledge of the potential created by in the 
far field? This is a classical inverse potential problem having in the general 
case an infinity of sohitions. One can imagine for example a charged ball 
with initially a spherically symmetrical but inhomogeneous charge density 
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0(0, r). It is clear that all the c„'s except Cq are zero in this case while (f){t, r) 
can be distributed arbitrarily. So we have an infinite number of solutions 
for satisfying the infinite number of constraints c„ = Cq 5o,n where 5o,n is 
the Kronecker-symbol, and so the same constants c„'s can correspond to 
different fields 4>{t,x) (see ref.[6]). 

Stefan problem as a singular limit of the "asymmetric" phase-field dy- 
namics: Now we consider the process described by the following equations: 

j^, + pUt = l V[((D+ + D_) + (D+ - VU] . . 

This equations encompass several well-known problems. If = D_ the 
system (23) coincides with the phase-field model [3] which describes a broad 
class of pattern formation problems where the order parameter (p is called a 
"phase" . and stand for the diffusions coefficients in the bulk where 
0=1 and = — 1 respectively (i.e. far from the "interface" defined as the 
surface described by 0(t, x) = 0). The first equation in the system (23) is a 
particular case of Eq. (1) for /9 = 0, and the second one is the time-dependent 
Ginzburg-Landau equation [7] in the presence of external field U. Here a 
and P are kinetic coefficients or the relaxation times of the fields and U 
respectively, ^ is the coherence length, and is the double-well potential 
depth. One can easily see that if a = (3 = and — D_ the system (23) 
is equivalent to the Cahn-Hilliard equation [8] 

0*= -V^(e'V^0+(0-0=^)/a) (24) 

describing the process of phase separation. It corresponds to the dynamics 
described by Eq.(9) when L = and the functional H has the Ginzburg- 
Landau free energy form: 

f .C^(V0)2 0^ 0\ ,n 

Jm^ 2 2a Aa' ^ ' 

As it was shown by Caginalp [4] the phase-field model (i.e. the Eqs.(23) 
with — DS) in the singular hmit ^ — > and a — > reduces to the 
symmetrical (i.e.when Dj^ — Stefan problem 

(3Ut = V^U xeM (26a) 

U = a{K + aVn) X e dQ{t) (266) 

Vn^-[dnU] xedn{t) (26c) 

where dD,{t) is a moving interface between the growing bubble fl{t) and the 
rest of the space, a is the surface tension proportional to i/^/a, k is the 
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mean curvature of the moving boundary dfl, Vn is the normal component of 
the local velocity of dQ{t), and [dn U] is the jump of the normal component 
of the gradient of U, dnU, across the boundary dO,{t). Laplacian growth 
corresponds to the quasistatic Stefan problem (or Hele-Shaw problem) when 
/3 = 0. 

The advantage of the "asymmetric" phase-field model (23) is that the 
same singular limit {C, ^ 0, a ^ 0) reduces the system (23) to the general 
(not necessarily symmetric) Stefan problem: 

/4 Ut = V^U xen± (27a) 

U = a{K + avn) X e dn (276) 

Vn = -D+dn U\phasel + D_dn U\phase2 X G dVL{t) (27c) 

Here /3± = I3/D±, while Q± are domains corresponding to phases (f) — ±1. 

One can further extend the system (23) to the more general processes 
described by equations 



<j)t + l3Ut = (pVqV + r)U 
ai^4)t^ ^^V^(l) + F'{(l))/a + U 



(28) 



Here p, q, and r can be functions of or of the space x. In other words 
the RHS of the first equation in the (28) is the arbitrary elliptic operator 
of the second order acting on U. In the RHS of the second equation of 
(28) the function -F(0) is introduced such that its derivative with respect 
to the scalar field has two stable solutions, 0±, and the unstable one, (po'- 
F'{4'±) = F'{4'o) = 0. It is easy to see that the same singular limit as used 
above reduces the system (28) to the following interface dynamics: 

P± Ut = (p±Vg±V + r±) U, x e Q± (29a) 

U^a{K + aVn) xedQ{t) (296) 

v^ = -p[q dn U] x e dQ{t) (29c) 

where p±, q±, and r± are the functions p, q, and r respectively defined 
in the domains and Q_. It was shown [9] that this kind of interface 
dynamics (describing rather "elliptic" than Laplacian growth, because of 
pVgV + r) possesses the infinite number of conservation laws for an 
arbitrarily dimensional process. 

Now we investigate the quantities c„ defined by Eq.(3) for the processes 
described by (28) in this singular hmit: 
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Conserved integrals Cn 's in the singular limit: The limit ^ — >^ makes 
the width between different phases infinitesimally thin, and the hmit a — > 
excludes fluctuations of (f){x), so it can take only three values: two stable, 
(f)±, and one unstable, 0o- Let us choose for simplicity F in the Ginzburg- 
Landau form: F' = (j) — (f)^. In this case (f)± = ±1, 0o — 0. The equation 
(j)Q{t,x) = describes the boundary moving between two phases = ±1. 
By introducing the new function 



we describe the interface x) = 1/2 moving between two phases described 
by 0' = 1 in Q{t) and 0' = outside it. Consequently, the integrals c„ in 
this limit are: 



But these integrals are the constants of motion found earlier in arbitrarily 
dimensional boundary motion when a = (3 = a = [9]. Thus, in the 
singular limit ^ — > 0, a — > the conserved quantities, c„, coincide with the 
conservation laws previously found. Therefore the field dynamics ((1) is a 
natural extension of the boundary motion d^{t) that possesses an infinite 
number of conservation laws. 

Integrability in the singular limit: An infinite number of polynomial con- 
servation laws is a necessary (but not always sufficient) condition for a 
nonlinear field problem to be exactly integrable [10]. When do the c„'s 
correspond to an integrable case? While the knowledge of the c„'s is not 
enough to recover the field 0(a;), the situation is much better in the singu- 
lar limit described above where the field dynamics (f){t) is reduced to the 
boundary dynamics dQ,{t). In this case one might recover the interface and 
the problem might be solvable. Indeed here the "matter distribution" 0' is 
especially trivial: cf)' — 1 in n{t) and (j)' — outside it. Therefore looking 
for the (j){t, x) by knowing the far-field potential $ is reduced to the looking 
for the boundary dQ(t) between two phases. This problem is called the 
"inverse potential problem" [11,12] and sometimes can be explicitly solved. 
For D = 2 and L = using conformal mapping it is possible to recover 
the boundary dQ from the knowledge of the singularities of the function 
Y^'^=i CnZ~'^ analytic in B? /VL [11], and where c„ are integrals defined by (3). 
However, if D 7^ 2 or L 7^ we do not have a constructive way to recover 
the boundary 9Q, but if Q, is (a) single-connected and (b) star-like (i.e. a 
point inside Vt exists such that all straight rays from this point cross the 
boundary dil only once) there exists a unique domain Q for given potential 



(/)' = (0 + l)/2 



(30) 




(31) 
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[12]. Nevertheless it is still unclear how to parametrize the corresponding 
boundary in this case. So, although the knowledge of c„'s in principle is 
enough to recover the boundary of a single-connected star-like domain in the 
general multi-dimensional case, only the 2-D Laplacian case is integrable in 
a constructive way (and not even for all kinds of singularities of the function 

Integrals Cn 's in the case of non-zero surface tension cr = C,/ \fo and non- 
zero kinetic coefficient a: If in the singular limit ,^ ^ 0, a ^ the ratio 
C/y/a has a finite non-zero value that corresponds to a non-zero surface 
tension, the integrals c„ defined as in Eq.(3) generally speaking are not 
conserved for the manifold M — i?^.(They are even not necessarily finite.) 
For the process described by the system(23) with P — the rate of Cm 
defined by (3) is 

d Cm/d t = D_ [ _ {ijjmdnU-Udn ll^m) dS""-' (32) 

where S^~^ is the hypersphere in whose radius tends to infinity when 
ipm diverges at infinity and tends to zero when ipm diverges at the origin. 
(Without surface tension it corresponds to the "internal" and "external" 
onc-sidc Stefan problems respectively, i.e. when U governed by the equation 
LU = is non-trivial only inside or outside the boundary dQ, while in the 
complementary phase U vanishes). If L = the function ipm diverges (at 
infinity or at the origin) as an m*'* power of r (see Eqs.(4)-(6)). In this case, 
for the RHS of Eq.(32) to vanish, U must decay where ipm diverges (i.e. 
at infinity or at the origin) faster than any power of r. Without surface 
tension, the field U which is zero in the phase where ipm diverges, since U 
is a harmonic function zero at the boundary, U — 0. Thus in this case the 
RHS of Eq.(32) is zero, so that the c^'s are conserved. 

For non-zero surface tension this is not the case, because U, satisfying 
non-trivial boundary conditions (see Eq.(27b)) decays at infinity (or at the 
origin) only algebraically, instead of exponentially (generally speaking as r 
at the origin or as at the infinity). Therefore, in this case, the RHS 
of Eq.(32) is not zero, so the c„'s are not conserved. For the 2D Laplacian 
internal problem (when ipn — z'") it is 

d Cn/d t = -4:D_TTnUn (33) 

where 

oo 

U{r, 6) - ^(r-"(C/„exp(-me) + exp(me))) (34) 

n=l 
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[13]. By the same reasoning the c„'s are not conserved when the kinetic 
coefficient a is not zero. 

However, there are situations when dcn/dt — even for non-zero sur- 
face tension and kinetic coefficient a: (i) First, this is the case if instead of 
L = the elliptic operator L is such that its solutions decay exponentially 
at infinity (for example if L = — a with a positive a). Evidently the 
surface integral in the RHS of Eq.(31) also vanishes, (ii) Second, 
conserved for the one-side (i.e. when = 0) non-zero surface tension case 
(which is beyond the phase-field model dealing only with the symmetrical 
two-side model, when = The one-side problem with non-zero sur- 
face tension is very common for different physical systems, for example for 
bubbling in Hele-Shaw cell, in solidification when one can neglect diffusion 
in a solid, etc. [2]. (For simplicity we consider here only isotropic surface 
tension. For the anisotropic case all the results hold, but a can no longer 
be pulled out of integrals and derivatives since it is no longer a number but 
instead an angle-dependent function.) 

It is easy to obtain the rate of change of c„ in this case directly from the 
boundary motion's formulation (29). Indeed because the integrand in (3) 
does not depend on time, but only on the domain of integration, we have 

dcjdt^ l/dt{[ - I )^mdV 

Jdn(t+dt) Jdn{t) 

(here dV stands for the element of the moving surface dfl{t). Further, it 
equals 

= / VnlpmdV = IprndnUdT 

Jdn(t) Jdn(t) 

because of Eq.(29c) (here p — 1, g+ = 1, and q- — are chosen for the 
simplicity) ; 

= / ijjrndnUdr - / {U - a {k + aVn)) Onlpm dV 

Jdn{t) Jdciit) 
because of Eq.(29b), and finally 

^ -cr {k+ adnU)dnipmdr (36) 
Jan 

since the vector field {U V ■0m — V'm V U) is divergenceless because of 
Eq.(29a): 

V(UV^|Jm - ^mVU) = 

So this is the law of change of Cm defined by (3) for a one-side Stefan 
model in the presence of surface tension and kinetics on the moving bound- 
ary. It is known [14] that without sources of the field U (i.e. when p — 0) 
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the boundary dfl{t) usually has the tendency to take a circular shape in the 
long-time asymptotics (excluding possible topological changes which are not 
considered here). In this case we expect d\c'^\/dt to he negative since |c^| 
indicates a deviation from a circular shape. It would be interesting to check 
this statement in the near future. 

Conclusions: We have shown straightforwardly that there exist an in- 
finite number of conserved quantities for the field dynamics described by 
Eq.(l), such that in the case of a bistable scalar field the singular limit 
corresponding to the boundary motion between two phases (Ising limit in 
the kinetics of phase transitions) can be integrable (for instance ii D — 2 
and L — V^). The criterion of integrability is the solvabihty of the as- 
sociated inverse potential problem, where the conserved integrals c„'s are 
the coefficients of the multipole expansion of the fictitious potential created 
by the matter distribution (f)(x,t). In the presence of surface tension and 
boundary kinetics these conservation laws are broken, but in the framework 
of the one-side model, the evolution of quantities c„'s has been calculated 
(Eq.(35)). It means that in many cases when the integrals c^'s completely 
determine the moving cluster, we can follow the interface' evolution with 
non-zero surface tension. 

I am grateful to F. J. Alexander, R. Almgren, and I. M. Gelfand for very 
helpful discussions of this problem and to S. Ponce-Dawson and P. Sievert 
for their comments on the paper. 
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